Asymptotic normalization coefficients from ab initio calculations 
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We present calculations of asymptotic normalization coefficients (ANCs) for one-nucleon removals 
from nuclear states of mass numbers 3 < A < 9. Our ANCs were computed from variational Monte 
Carlo solutions to the many-body Schrodinger equation with the combined Argonne wis two-nucleon 
and Urbana IX three-nucleon potentials. Instead of computing explicit overlap integrals, we applied 
a Green's function method that is insensitive to the difficulties of constructing and Monte Carlo 
sampling the long-range tails of the variational wave functions. This method also allows computation 
of the ANC at the physical separation energy, even when it differs from the separation energy for 
the Hamiltonian. We compare our results, which for most nuclei are the first ab initio calculations of 
ANCs, with existing experimental and theoretical results and discuss further possible applications 
of the technique. 
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Substantial experimental and theoretical effort over 
the past decade and a half has been expended on the ex- 
traction of asymptotic normalization coefficients (ANCs) 
from experiments involving light nuclei [ll-[9|. Most of 
this work has been motivated by the connection between 
ANCs and astrophysical cross sections, but ANCs also 
offer opportunities for significant new tests of ab initio 
nuclear calculations. In this Rapid Communication, we 
present predicted ANCs for several states of light nuclei 
up to ^ = 9, using the variational Monte Carlo (VMC) 
method and a realistic Hamiltonian. 

Recent years have seen rapid advances in the ab initio 

ut- 



theory of light nuclei 



ing power has been brought to bear on the problem of 
computing properties of light {A < 12) nuclei from a 
new generation of accurate nucleon-nucleon and three- 
nucleon potentials. Many nuclear properties have been 
computed from the modern nuclear interactions, includ- 
ing charge radii, electroweak transition amplitudes, cross 
sections for scattering and radiative capture, and spectro- 
scopic factors. Some ANCs have been computed 13-13], 
but there has been no systematic ab initio investigation 
of ANCs. 

An ANC characterizes the asymptotic form of a nuclear 
overlap function, which is the projection of a nuclear wave 
function onto a product of subclusters. We consider only 
cases of one-nucleon removal, so the subclusters within a 
nucleus of mass A are the removed or "last" nucleon itself 
and a residual nucleus of mass A — 1. (Although we refer 
to the "last nucleon," our wave functions are explicitly 
antisymmetric.) The overlap channel is further specified 
by orbital angular momentum / and its vector sum j with 
the spin of the last nucleon. The overlap function is then 



Ja 



(1) 

where ^'^^ is the wave function of the mass-A nucleus 
with angular momentum J a, '^a^J^i ^ specific state 



of the residual nucleus with angular momentum Ja-i, 
X is the spin-isospin vector of the last nucleon, and r^c 
is its separation from the center of mass of the other 
A — I nucleons. Square brackets denote angular momen- 
tum coupling, Yi are spherical harmonics, and A antisym- 
metrizes the product '^'^SiX^i with respect to particle 
exchange. The integral extends over all particle coordi- 
nates R = (ri, r2, r^). 

The form of the overlap as r — >■ cxd is well known, be- 
cause it satisfies a one-body Schrodinger equation includ- 
ing at most a Coulomb interaction. This form contains 
a Whittaker function W-rim- 



Ri^^-''Hr 



oo 



(2) 



with rj = aZA-iZN\/ fic^ /2B, a the fine-structure con- 
stant, Za~i and Zn respectively the charges of the resid- 
ual nucleus and the last nucleon, fi their reduced mass, B 
the separation energy of the last nucleon, k = y/2iiB/h, 
and m — I + 1/2. (If the last nucleon is a neutron, then 
Zn = and W-Tfrn{2kr) — ^2kr/TrKm{kr), a modified 
spherical Bessel function of the third kind.) In the fol- 
lowing, we omit the labels J a and Ja-i for compactness 
of notation. 

The only quantity in Eq. ([2]) that is not determined 
fully by the quantum numbers and the corresponding 
separation energy is the constant Cij. It characterizes 
the overall scale of the long-range A-body wave function 
in the Ij channel, and it is the ANC of that channel. 

It can be shown that although the spectroscopic factor 
Sij = J Rfj{r)r^dr may depend strongly on the short- 
range potential and the choice of wave-function repre- 
sentation, the ANC as both a theoretically-computed and 
an experimentally-inferred quantity is less dependent on 
such details [ll,[3|- Given reactions (e.g., well below the 
Coulomb barrier) that probe only the asymptotic part of 
Rij, ANCs can be extracted from data with fewer as- 
sumptions than spectroscopic factors can. 

Computing an ANC by direct integration of Eq. ([T]) 
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FIG. 1: (Color online) Points with Monte Carlo statistical er- 
rors show the *Li n ^Li overlap, computed from Eq. Q, of 
our VMC wave functions in pi/2 (red squares) and P3/2 (blue 
circles) channels. Curves with error bands show the asymp- 
totic forms in Eq. scaled by ANCs from Eq. ([5]). Dashed 
(red) curves are asymptotics for pi/2 and solid (blue) ones 
are asymptotics forp3/2- They are labeled "Bh" and "Bcxpt" 
according to the assumed neutron separation energies. 



is problematic for most many-body methods. First, ab 
initio calculations may not yield the correct asymptotic 
form of Eq. ([2]). For example, methods using a harmonic- 
oscillator basis have basis functions with an asymptotic 
form e"'-''/^-' , so that convergence to a long-range asymp- 
totic form similar to e"'^'" is slow. In variational methods, 
it is often difficult to construct a consistent set of correla- 
tions that has good long-range asymptotics while retain- 
ing short-range properties that are important for the vari- 
ational energy. Second, the assumed Hamiltonian may 
not reproduce the experimental separation energy -Boxpt 
even when wave functions arc computed exactly. Third, 
Monte Carlo methods suffer from the difficulty of find- 
ing a sampling scheme that samples the tails of Eq. ([T]) 
thoroughly while minimizing sample variance. All three 
difficulties are illustrated in Fig. [1] 

There is another approach to computing ANCs that 
avoids all three of these difficulties, and versions of it 

23-Ei- In this 



have been derived in several contexts 
approach, explicit computation of the overlap function is 
replaced by an integral over the wave-function interior. 
The Schrodinger equation 



{H-E)^A = 



(3) 



that yields wave function a with energy E may be 
rewritten as 

= -[T,,^ + Vc + B]-\u,,^-Vc)^A 

- [T,el + Vc+B]-^ (i7i„t - i^int ) ■ (4) 

We have broken up the Hamiltonian H into the rela- 
tive kinetic energy Trci between the residual nucleus and 



last nucleon, a sum of terms i?int involving only nu- 
cleons within the residual nucleus, and a sum of terms 
Ural involving the last nucleon. The point-Coulomb po- 
tential between the residual nucleus and last nucleon is 
Vc = ZA-iZ]^ahc/rcc- Similarly, E = E^t — B, with 
Eint being the purely internal energy of the residual nu- 
cleus. 

If we rewrite the Green's function [T^-d + Vc + B]^^ 
in terms of special functions, project onto the product 



A-i [xYim^ 
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identity that {Hint ' 
limit, we find that 
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as in Eq. ([T]), take advantage of the 
^^int)^'A-i = 0, and take the r -> 00 

(5) 

-^'LlX^^'/Cfcc) (f/rel - Vc) ^AdK. 



The integral extends over all particle coordinates, M^rim 
is the Whittaker function that is irregular at infinity, 
w is its Wronskian with the regular Whittaker function 



and the angular momentum algebra is omitted 



for simplicity. 

The utility of Eq. 



arises from the form of Urc\- If 



Vij and Vijk are respectively terms of the two- and three- 
body potentials involving nucleons labeled i,j, and fc, 
and we always label the last nucleon A, then 



rcl 
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(6) 



At large separation r^c of the last nucleon, only the 
Coulomb terms of ViA are nonzero. The monopole term of 
their sum is equal to Vc: so the factor U^c\ — Vc in Eq- ® 
is short-ranged. In our calculations, it limits significant 
contributions to r^c < 7 fm. Equation ([5]) thus reduces 
a problematic calculation involving the outer regions of 
to a manageable calculation involving its interior. 

We implemented Eq. ([5]) within the VMC method de- 
scribed in Ref. 2J]. The Hamiltonian comprised the Ar- 
gonne vis two-nucleon [25| and Urbana IX three-nucleon 
interactions 26|. For this interaction (AV18-I-UIX) we 
constructed variational wave functions and 
that minimize the energy expectation values while con- 
straining them to give approximately correct charge radii, 
as determined experimentally (where known) or by ex- 
act Green's function Monte Carlo (GFMC) calculations. 
The ANC integral was performed by Monte Carlo inte- 
gration, using the same sampling scheme (with weight 
proportional to I^Ap) as our energy calculations. 

The distribution in Vcc of the ANC integrand is shown 



in Fig. [2] for the specific case of *Li 



^Li. (Where 



there is no further label, the ground state of a nucleus is 
implied). It may be seen that the ANC integral is con- 
tained entirely within about 7 fm. The distribution of 
Monte Carlo samples, shown as a dotted curve, is broadly 
similar to the distribution of the ANC integrand, so the 
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FIG. 2: (Color online) The integrand of Eq. © {x2^/khw) 
is shown for the pi/2 (red squares) and p-j,/2 (blue circles) 
neutrons in *Li n^Li. It is binned by the n-^Li separation 
Tec with bars showing Monte Carlo errors. The solid curves 
are cumulative integrals of Eq. ^ , starting from the origin; at 
large Vcc, they are the ANCs (divided by 2 for visibility on this 
scale) . The dotted curve with no scale shows the distribution 
of Monte Carlo samples. 

integral is computed with relatively small statistical er- 
rors. 

The computed Cij depend sensitively on the separation 
energies B. Equation ([5|) contains B implicitly through 
k = \/2jlB/h and 77 cx 1/^fB, and it is rigorously true 
when B = E-mt — E for the given potential. However, 
there are often significant differences between this B and 
the experimental separation energy iJexpt • We computed 
several ANCs in light nuclei, first using the GFMC Bh 
for the AVIS+UIX Hamiltonian and then using Boxpt- 

The use of i?Gxpt in Eq. ([5]) may be understood by con- 
sidering small changes to the potential. When B <^ \E\, 
they can produce small changes in the wave-function 
interior but large fractional changes in B. The short- 
range part of the variational wave function derived from 
AV18-I-UIX is, therefore, similar to the solution that 
would be obtained from a slightly different potential 
tuned (e.g. with small extra terms) so that Bh = Sexpt- 
Inserting a fc oc ^i^oxpt hito Eq. ([5]) matches a wave- 
function interior approximating the true wave function 
onto the asymptotic form corresponding to -Bcxpt. In- 
structive illustrations of this general principle, applied to 
much simpler wave functions, may be found in Ref. [27| . 

The use of Eq. ([5]) to compute asymptotic overlaps 
is demonstrated in Fig. [l] where ^Li — ^ n^Li over- 
laps computed directly from Eq. ([T]) are plotted next to 
CijW^r^m/r from Eq. ([SJ. It can be seen that the W^nm 
corresponding to Bh = 1.3 MeV are rather different 
from those for Sexpt = 2.03 MeV, though both energies 
are small fractions of the 41.3 MeV total binding energy 
for ^Li. 

For both B values, the asymptotic Rij match the short- 
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FIG. 3: (Color online) Predicted ANCs from Eq. (O, divided 
by experimentally-derived values from the references given at 
the right (those not appearing elsewhere are Refs. [29l - [3ll |). 
For each ANC, small error bars indicate the Monte Carlo error 
of Table U and lar ger error bars indicate its quadrature sum 
with the experimental error. Results for the same computed 
ANC divided by different "experimental" numbers are joined 
with dashed lines. Parentheses indicate particularly uncertain 
experimental constraints. 



range overlaps at ^ 4 fni, where the ANC integral starts 
to converge. Use of Bh yields C^^/j = 0.029(2) fm"^ 
and Cp3/2 = 0.237(9) fm~^, compared with the respec- 
tive values 0.048(6) fm~^ and 0.384(38) fm"^ from a 
transfer- reaction study . The match between the com- 
puted and "measured" results is poor. Using Bexpt yields 
0.048(3) fm"^ and 0.382(14) fm-\ in very good agree- 
ment with experiment. This pattern of agreement with 
experiment for i?cxpt but disagreement for Bh repeats 
in all cases of substantial difference between Bh and 
Bcxpt- In the following, we consider only ANCs computed 
from -Bcxpt, and we assign uncertainties based entirely on 
Monte Carlo statistics rather than (difficult) assessments 
of the variational wave functions. Limited testing with 
variant wave functions suggests that the total uncertainty 
is not much larger than the statistical uncertainties. 

Our results are shown in Table [J and compared 
with experimentally-derived numbers (where available) 
in Fig. [21 The lowest three sections of Table |T] re- 



4 



TABLE I: ANCs computed from Eq. ((5| for given yl-body 
nuclei, {A — l)-body residual nuclei, and angular momentum 
channels Ij or Units are fm-i/2^ and / -wave ANCs 

have been multiplied by 10"^. Error estimates reflect Monte 
Carlo statistics only, and columns left empty are zero by exact 
symmetries. Asterisks denote first excited states. 



A 


A - 1 


Sl/2 


d3/2 


Cd 3/2 


/Cs 1/2 


■^H 




2.127(8) 


-0.0979(9) 


-0.0460(5) 


^Hc 




2.144(8) 


-0.0927(10) 


-0.0432(5) 


■'He 




-6.55(2) 








■*He 


^He 


6.42(2) 








A 


A-1 


Pl/2 


P3/2 


h/2 X 10^ 


/7/2 X 10^ 


''Li 


«He 




3.68(5) 






''Li* 


''He 


3.49(5) 








^Li 


''Li 


1.652(12) 


1.890(13) 


-78(20) 




^Li* 


"Li 


-0.543(16) 


-2.54(4) 






^Bc 


"Li 


-1.87(3) 


-2.15(3) 


63(9) 




^Bc* 


"Li 


0.559(16) 


2.59(5) 






**Li 


"Li 


0.218(6) 


-0.618(11) 


5.2(5) 


-2.5(15) 


**Li* 


"Li 


-0.090(3) 


0.281(5) 


- 0.6(2) 




**B 


"Be 


0.246(9) 


-0.691(17) 


1.1(2) 


-1.1(5) 




**B 


-0.309(7) 


1.125(12) 


1.9(5) 


-0.5(18) 




**Li 


0.308(7) 


-1.140(13) 


- 4.1(10) 


5(3) 


''Li 


**Li* 


-0.122(3) 


0.695(7) 


-1.1(6) 




9t ■ 

Li 


8tt 

He 




— 5.99(8) 






''Bc 


**Li 


5.03(6) 


9.50(11) 


35(34) 


257(112) 


Be 


'^Li* 


6.56(5) 


-6.21(7) 


364(40) 




A 


A — 1 


'i 

P 


4 

p 


'^f X 10^ 


"/ X lO-' 


''Li 


"Li 


2.510(18) 


0.029(18) 




- 78(20) 


"Li* 


"Li 


-2.57(5) 


-0.33(3) 






"Be 


«Li 


-2.85(4) 


-0.04(4) 




- 63(9) 


"Be- 


''Li 


2.63(5) 


0.34(3) 






^Li 


8Li- 


-0.599(7) 


-0.373(7) 




1.1(6) 


^Be 


«Li* 


-0.25(9) 


-9.03(8) 




-364(40) 


A 


A-1 




'p 


*/ X 10^ 


■^Z X 10-' 


'■>c 


«B 


0.868(14) 


0.779(12) 


0.1(19) 


-2(1) 


^Li 


«Li 


-0.882(15) 


-0.785(12) 


3.3(34) 


5.2(19) 


^Be 


'^Li 


10.75(12) 


-0.25(10) 


256(117) 


42(65) 


A 


A-1 


''p 


'p 


■'f X lO'' 


■'f X 10-^ 


«Li 


'Li 


-0.283(12) 


-0.591(12) 


- 0.3(16) 


- 5.8(10) 


»Li* 


"Li 


0.220(6) 


0.197(5) 




0.6(2) 


»B 


"Be 


-0.315(19) 


-0.662(19) 


- 0.6(5) 


- 1.4(4) 



peat information from the second section, but in "chan- 
nel spin" coupling of the form /] , instead 



Ja-1 



Ja 



We examined most channels up to 



of 

A = 9 with either the ^-body or the residual nucleus 
in its ground state and with both stable against particle 
decay. We now comment briefly on the comparison of 
our results with past work. Extensive discussions of past 
experimental and theoretical estimates may be found in 
Refs. [13,132 1 . 



The s-wave ANCs for A < 4 nuclei have typically been 
inferred from cross sections using techniques based on 
analyticity of the scattering amplitude [H, [s^l, mostly 
thirty or more years ago. Although our ANCs agree with 
many of those results, Fig. [3] demonstrates the consider- 
able systematic uncertainties of those methods discussed 
in Refs. 



3al34| 



ANCs of '^H and '^He have been computed previously 
from modern realistic interactions using Eq. ^ [l^ and 
were the focus of much activity following the development 
of Faddeev methods (s^ - fsij . Particular emphasis was 



placed on the ratio Cd3/2/Cs i/2jinost precisely inferred 
from tensor analyzing powers [36|; those results are in 
reasonable agreement with ours. 



The Pisa group has computed ANCs for ^ < 4 |15l.ll6| 
with AV18-FUIX. Their A = 3 Csi/2 are within 0.5% of 
ours, but their 0^3/2/0^1/2 are 10% smaller. Their "^He 
ANCs are also about 6% smaller than ours. The reason 
for this difference is unclear; it could reflect shortcom- 
ings of the variational wave functions, which miss the 
true AV18-MJIX binding energy by 850 keV in ''He. On- 
going work to compute overlaps using essentially exact 
wave functions from the GFMC method seems to sup- 
port our values of the A = A ANCs 37|. (Nuclei with 
A = 3, 4 have substantially identical ANCs for Bh and 
Scxpt because the AV18-I-UIX interaction was tuned to 
have Bh — i?expt in these systems. Pisa ANCs converted 
to our conventions may be found in Ref. [22].) 

For A > A ANCs, experimental constraints have been 
inferred almost entirely from transfer [iKl, [3, [111, 
knockout or breakup |6l| reactions, and are of gener- 
ally more recent vintage than the A < 4 ANCs. In some 
cases components of different j contribute indistinguish- 
ably to differential cross sections, which then constrain 
only the sum J^j Ofj- These cases are indicated in Fig. [3] 
and shown as the square root of the sum for comparabil- 
ity of error bars. Our p-shell ANCs are in broadly good 
agreement with those inferred from experiment, particu- 
larly for the well-measured A = 8 ground state ANCs as 
discussed above. (Our calculations for A = 8 also agree 
with prior theoretical estimates of [13, [i^-) Reference 



27j presented many ANCs computed by applying Eq. 
([5]) with a simpler potential to harmonic-oscillator wave 
functions derived from shell models; about half of our 
p-shell ANCs disagree with those calculations by more 
than 25%. 

The most significant differences from previous work are 
in the ''Li n ^Li ANCs. The comparison with exper- 
iment here is difficult because of the wide range of esti- 
mates, which extend from Cf^ = 1.26 to 2.82 fm~^/^ 

just from {d,t) at varying energy ([3], with full range 
shown in Fig . [3]) and include other values within that 
range [H, [io]. The effective ANC of Huang et al. [ill , 
whose capture model successfully matches ^Li{p, 7)^Be 
data, is 25% below ours. 

The theoretical ANCs for ^Li — >■ n ^Li (from a simpler 
model) in ReL (27)] are 20% to 40% smaller than ours. As 
with ^He, ongoing GFMC work (with an improved three- 
body interaction) seems to support our results [st']. We 
also disagree with earlier integral-method predictions of 
the ratio of ''Be — > p^Li to isospin-mirror ^Li n^Li 
ANCs fi^l, finding 1.15 instead of 1.05 (though we agree 
with their 1.12 as the ratio of to ^Li ANCs). The 
sources of these differences are unclear. 

Table [J includes ANCs for both p- and /-wave chan- 
nels of p-shell nuclei. The small /-wave components arise 
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from the tensor terms the Hamiltonian, analogously to 
the cZ-wave components in s-shell nuclei [43j . We are un- 
aware of any previous calculations of /-wave ANCs or at- 
tempts to measure them. A DWBA calculation of tensor 
analyzing powers in sub-Coulomb ^'^^'^°^Pb(^Li,^ Li)X 
(analogous to triton d/s ratio experiments) suggests 
that both cross sections and analyzing powers may be 
too small to allow measurement of the f/p ratio [28| . 
Nonetheless, the /-wave ANCs demonstrate something of 
the power of the integral method: Within the VMC ap- 
proach, computing ANCs for these small- amplitude chan- 
nels from Eq. ([1]) would require far more computing time 
to achieve the same statistical accuracy, even if our varia- 
tional wave functions guaranteed the correct asymptotic 
form. 

Several extensions of this technique within the context 
of quantum Monte Carlo methods suggest themselves. 
The overlaps need not correspond only to one-nucleon 
removal, but may include cluster overlaps like *He — ^ dd 
and ^Be — ^ a'^He. A straightforward extension of the 
definition of ANCs to include unbound states allows the 
prediction of energy widths from the integral method [S- 
|. The ANC integral can also be evaluated within the 
GFMC method, which provides essentially exact results 
for a given potential. Use of the (computationally more 
demanding) Illinois three-body potentials ^] to generate 
the wave function and/or the ANC kernel will provide 
more accurate ANCs and Bh closer to i3oxpt- Finally, 
use of Eq. Q away from the r ^ oo limit should allow 
more accurate calculations of overlaps at all radii [13, [SH, 

13, m. 
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